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Abstract Critical Rayleigh numbers have been measured in 
a liquid metal cylinder of finite height in the presence of a 
rotating magnetic field. Several different stability regimes were 
observed, which were determined by the values of the Rayleigh 
and Hartmann numbers. For weak rotating magnetic fields and 
small Rayleigh numbers, the experimental observations can 
be explained by the existence of a single non-axisymmetric 
meridional roll rotating around the cylinder, driven by the 
azimuthal component of the magnetic field. The measured 
dependence of rotational velocity on magnetic field strength is 
consistent with the existence of laminar flow in this regime. 

List of symbols 

B 0 magnitude of magnetic induction 

B n B 0 radial and azimuthal magnetic induction com- 

ponents 

C wall admittance 

d cell diameter 

d w wall thickness 

g gravity at earth’s surface 

Ha Hartmann number 

h cell height 

kj thermal conductivity of fluid 

k w thermal conductivity of wail 

LI, L2, L3, L4 thermistor temperatures 
Ra Rayleigh number 

Ra‘ critical Rayleigh number for the transition from 

no flow to laminar flow 

Ra* critical Rayleigh number for the transition from 

time-independent to time-dependent flow 
r radial coordinate 

T a temperature at top of cell 

T b temperature at bottom of cell 
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AT temperature difference between cell bottom 

and cell top 

AT C critical temperature difference between cell 

bottom and top time 
t time 

U \ , (72, (73, 

(74 thermistor temperatures 

z vertical coordinate 

a volumetric thermal expansion coefficient 

8 skin depth 

k thermal diffusivity 

ft magnetic permeability 

v kinematic viscosity 

p density 

o electrical conductivity 

0 azimuthal coordinate 

a) angular frequency of magnetic induction 

1 

Introduction 

There has been considerable interest lately in the application 
of rotating magnetic fields to the processes of single crystal 
growth (see, for example, Dold et al. 1995; Bruckner and 
Schwerdtfeger 1994; Barmin et al. 1992). This has been 
motivated by the potential benefits that rotating magnetic 
fields can have on the crystallization process. These 
include stirring of multicomponent melts, enhancement of the 
interface stability, and a reduction of thermal asymmetry 
at the solidification front. In order to take advantage of 
these benefits, it is essential to understand the fluid flow 
of liquid metal under the influence of a rotating magnetic 
field. Indeed, there have been a number of papers devoted to 
this endeavor (Gelfgat et al. 1993; Davidson and Hunt 1987, 
Doronin et al. 1973). This paper presents experimental results 
of fluid flow behavior under the combined forces generated 
by rotating magnetic fields and gravity. To achieve this, 
a cylindrical test cell has been designed to make carefully 
controlled measurements of the flow transitions between the 
regimes that exist under these combined forces. These flow 
transitions are observed by measuring the temperature 
changes of thermistors immersed in the liquid metal. Critical 
Rayleigh numbers are measured for the case of B 0 = 0 
for comparison to previous theoretical, numerical, and 
experimental results. Then, the critical Rayleigh numbers 
are measured in the presence of a rotating magnetic field, 
and the fluid behavior in the flow regimes identified are 
examined. 
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Experimental set-up 

The experiments were conducted using liquid gallium con- 
tained in a vertical cylindrical test cell depicted in Fig. 1. The 
side walls were made from polymethyl methacrylate (PMMA) 
and copper discs were used for the top and bottom. The tem- 
peratures of the top and bottom of the cell were controlled by 
fluid recirculating from constant temperature water baths. The 
transmission of vibrations between the fluid pumps and the 
test cell were reduced by first flowing the fluid through passive 
pressurized containers placed between the fluid pumps and the 
test cell. The temperatures of the cell top ( T a ) and cell bottom 
( T b ) were measured by thermocouples placed in the copper 
disks less than 0.5 mm from the copper gallium interface. The 
cell height was 24 mm, and the cell diameter was 23 mm for 
a resultant aspect ratio hid of 1.04. Eight thermistors were 
placed in the cell to record temperature changes and infer from 
them the nature of the fluid flow. Four of the thermistors, 
labelled LI through L4, were placed in a horizontal plane at 
z — hl 3, and the other four, labelled U\ through U4, were placed 
in a horizontal plane located at z = 2hl3. The thermistors were 
positioned at azimuthal positions of 0 = 0, n!2, n, and 37t/2. The 
diameters of the thermistors were 0.6 mm and they protruded 
into the melt a distance of 2 mm. The thermistors were 
operated with a constant current of 10 p A which generated 
negligible self-heating in the system. Both the thermocouples 
and thermistors were calibrated by recording their voltages 
under isothermal conditions at several temperatures spanning 
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Fig. 1. Schematic drawing of the test cell. Shown at the top is a vertical 
cross section of the cell in the r, z plane for 0 = 0. The gravitational 
vector is in the -z direction. At the bottom is a horizontal cross 
section of the cell in the rO plane for z = 2fi/3. LI, L2, L3, L4 and 171, U2, 
(73, U4 designate the thermistors immersed in the liquid gallium 
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Fig. 2. A side and top view of the position of the test cell within the 
rotating magnetic field coil 


the temperature range of interest. The thermocouple and 
thermistor signals were digitized with a 16-Bit analog-to-digital 
convertor (ADC) located in a slot of an AT-type PC and then 
recorded on computer disk. Up to 8 thermistor or ther- 
mocouple signals could be recorded simultaneously. The noise 
of the system was such that it was possible to obtain a temp- 
erature resolution of better than 2 mK with the thermistors 
immersed in liquid gallium. 

The test cell was placed inside a rotating magnetic field as 
sketched in Fig. 2. The field was generated by a three pole, three 
phase inductor operating at 60 Hz. The magnet had a 120 mm 
inner diameter, a 180 mm outer diameter, and was 120 mm in 
height. The magnetic field profiles were measured with a Hall 
probe without the test cell present and the variation in field 
strength over both the vertical and horizontal dimensions of 
the cell was less than 3%. The magnetic field rotated in the r, 0 
plane and can be described by the equations: 

B r ( t ) = B 0 cos (0- cot), B 0 (t) = B o sin{O-cQt) ( 1 ) 

where B Q is the magnitude of the applied magnetic field and 
co is its angular frequency. 
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Measurement of critical Rayleigh numbers 

The test cell was designed to accurately measure critical values 
of the Rayleigh number. The Rayleigh number is defined as 

Ra= 8 a{Tb ~ T ‘‘ ) - ( 2 ) 

KV 

where g is the gravitational acceleration on the earth’s surface, 
9.81 m/s 2 , x the volumetric thermal expansion coefficient, 
T b — T a the vertical temperature difference across the cell, h the 
cell height, k the thermal diffusivity and v the kinematic 
viscosity. Table 1 lists the material properties of liquid gallium 




Property (323 K) 

Symbol 

Value 

Units 

Reference(s) 

Thermal expansion 

3t 

1.25x10 4 

1/K 

Hoather (1936) 

Thermal diffusivity 

K 

1.38 x 10~ 5 

m 2 /s 

Iida and Guthrie (1988) 
Powell and Childs (1974) 

Kinematic viscosity 

V 

3.177 x 10 7 

m 2 /s 

Grosse (1961) 

Density 

P 

6.075 x 10 3 

kg/m 3 

Nalgiev and Ibragimov 
(1974) 

Electrical conductivity 

(7 

3.82 x 10" 

S/m 

Cusack and Kendall (I960) 


at 323 K necessary to calculate the Rayleigh number and the 
interaction of liquid gallium with a rotating magnetic field. The 
thermal diffusivity was calculated from the thermal conductiv- 
ity given by Powell and Childs (1972) and the heat capacity 
found in Iida and Guthrie (1988). The material properties a, k, 
and v vary as a function of temperature. However, the quantity 
7.1k v varies by less than 3% over the range of temperatures at 
which the critical Rayleigh numbers were measured. Therefore, 
all calculations of the Rayleigh number assumed that the 
material properties were temperature-independent and utilized 
the values listed in Table 1. 

The critical Rayleigh number Ra‘ is generally defined to be 
the transition value between the static condition state and 
the time-independent laminar convection state. The critical 
Rayleigh number Ra 1 is the value where the transition between 
time-independent steady convection and time-dependent con- 
vection occurs. These critical values were measured here to 
compare to previous theoretical results and as a baseline for 
comparison of these transitions in the presence of a rotating 
magnetic field. 

Two methods were used to measure the critical Rayleigh 
numbers. One method consisted of starting with a vertical 
temperature difference AT less than the critical AT , keeping 
T a constant, and slowly raising the temperature of the cylinder 
bottom T b . The second method consisted of starting with a zlT 
larger than the critical AT, keeping T a constant, and slowly 
lowering T b . As pointed out by Muller and Neumann (1983), it 
is necessary to keep the heating rate quasistatic, so that the 
temperature profile in the fluid remains linear. This condition 
can be met as long as 

dr 2AT c k 

d7 < ~ir (3) 

where dT/dt is the heating rate and AT C is the temperature 
difference T b — T a when the critical Rayleigh number is 
reached. For a low Prandtl number fluid such as gallium, this 
criterion is easily met, and for the present experimental set-up, 
this criterion dictates that dT/dt <5 x 10 2 K/s. 

The experimental results suggest that Eq. (3) is a necessary 
but not sufficient condition to obtain an accurate measurement 
of Ra 1 . There is an additional condition which occurs as a result 
of fluid inertia. Consider a fluid at rest and placed in the 
temperature bath, so that Ra > Ra c . Then the rest state is 
unstable, and it will relax to the convection state. The time 
constant of this transition is approximately proportional to the 
difference Ra — Ra 1 . Retardation effects for the fluid flow 
rearrangement will be seen for any finite temperature ramping 



Fig. 3. Temperature differences of thermistors situated directly across 
from each other in the test cell versus Rayleigh number 


rates. For B 0 = 0, hysteresis was observed in the measurement 
of Ra c for heating and cooling rates \dTldt\ >2 x 10 ~ 3 K/s. That 
is, for heating rates larger than approximately 2 x 10 ~ 3 K/s, the 
onset of convection occurred at Rayleigh numbers above Ra ( , 
and for large cooling rates the transition to quiescence oc- 
curred below Ra For B 0 = 0, Experimental measurements of 
Ra c found it to be independent of heating rates below 
2 x 10 " 3 K/s. For B 0 ^ 0, hysteresis was observed in the 
measurement of Ra c , even for heating and cooling rates \dT/dt\ 
<4 x 10 “ 4 K/s. 

Figure 3 is a plot of the temperature differences of thermi- 
stors situated 180 across from each other versus Rayleigh 
number. For this particular experiment, the 8 signals recorded 
were from thermistors LI L4, C/1, £73, and thermocouples 
T a and T h . The break in the slope of L4 L2 indicates the onset 
of heat transport by convection and the occurrence of Ra c . It 
has been predicted (Buell and Catton 1983) that for a vertical 
cylinder with an aspect ratio of 1 the flow mode which will 
occur immediately above Ra L is a single non-axisymmetric 
meridional roll (a roll in the r-z plane). The data in Fig. 3 are 
consistent with this prediction. The data indicate that a single 
roll exists predominantly in the L2, L4 z plane. The C/3 — C/1 
and L3 — LI temperature differences remain nearly zero be- 
cause the lines drawn between these thermistors are perpen- 
dicular to the direction of flow. 

The measured value of Ra ‘ without a rotating magnetic field 
was 3800. As shown in Table 2, this compares favorably to 
recent numerical simulations (Neumann 1990) and theoretical 




Table 2. 


Ra 1 


Method 


Reference 


3800 

Experimental 

This work 

3610 

Numerical 

Neumann (1990) 

3800 

Theoretical 

Buell and Catton (1983) 

6100 

Experimental 

Weber (1988) 

6000 

Experimental 

Neumann (1986) 


results (Buell and Catton 1983). The numerical and theoretical 
results shown in Table 2 assume a vertical cylinder with an 
aspect ratio of 1 and adiabatic sidewalls. Catton (1972) 
demonstrated that the thermal conductivity of the swidewalls 
can have a significant effect on Ra\ He defined a wall ad- 
mittance C = hk f !d w k w y where k f and k w are the thermal 
conductivities of the fluid and the wall and d n , is the wall 
thickness. The wall admittaance of the present experiment 
(C ^ 330) is such that the experimental data can be compared 
to theoretical calculations in the limit of adiabatic sidewalls. 

Figure 4 is a plot of the derivative of the thermistor signals 
with respect to time versus Rayleigh number. The thermistors 
signals have been offset for clarity. By taking the derivative of 
the thermistor signals the transition to time-dependent convec- 
tion becomes quite pronounced. A hysteresis etfect was 
observed in the measurement of Ra l , even down to the lowest 
heating and cooling rates used, |d7Vdt| ^ 1.2 x 10 4 K/s. At the 
lowest cooling rate used, Ra { occurred at an approximately 
0.5 K lower temperature upon cooling than it did upon heating. 
Such hysteresis effects have been experimentally observed 
previously (see, for example, Krishnamurti 1973), and are most 
likely the result of finite amplitude instabilities. In the present 
experiment it was found that Ra ! ^ 7.3 Ra l , This is comparable 
to the experimental results of Miiller et al. (1984), who found 
Ra { ^ 8 Ra 1 '. However, it is unclear why the value of Ra ! upon 
heating, 2.7 x 10 4 , is significantly higher than the recent 
numerical result of 1.5 x 10 4 (Neumann 1990). 



Fig. 4. Derivatives of thermistor temperatures with respect to time 
versus Rayleigh number. The transition to time-dependent flow is 
observed at Ru = 2.7 x 10 4 
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Flow transitions in a rotating magnetic field 

Consider a cylinder of liquid metal such as depicted in Fig. 2 
under the influence of a magnetic field described by Eq. (1). 
The rotating magnetic field induces electric currents in the 
liquid which in turn interact with the field to produce an 
azimuthal electromagnetic force in the melt. This force drives 
an azimuthal fluid rotation in the cylinder. The strength of the 
magnetic field used in these experiments was such that the 
angular frequency of induced fluid rotation was always several 
orders of magnitude less than 60 Hz, the frequency of the 
applied magnetic field. The degree of penetration of the field 
into the liquid is determined by the skin depth h where 



where ft is the magnetic permeability and a is the electrical 
conductivity. Here, the skin depth is larger than the diameter 
of the test cell and so the transverse magnetic field can be 
satisfactorily approximated as uniform and independent of r. 

Two non-dimensional numbers are required to describe the 
fluid flow states that exist under the combination of buoyancy- 
driven convection and convection driven by a rotating 
magnetic field. These non-dimensional numbers can be chosen 
as the Rayleigh number, defined in Eq. (2), and the Hartmann 
number, defined as 


Ha = hB 0 



( 5 ) 


Experimentally, these numbers were varied by changing the 
temperature difference T b ~ T a and the magnetic field strength 
£ 0 , respectively. As in the case when £ 0 = 0, Ra c and Ra' are 
used to define the critical numbers at which the fluid flow 
transitions from one state to another. Howewver, an azimuthal 
component of fluid flow always exists in the presence of 
a rotating magnetic field. Therefore, Ra c is now defined as 
being that value at which buoyancy-driven meridional flow 
exists and contributes significantly to convective heat transfer 
in the cylinder. Ra ' is now defined as that value at which 
time-dependent flow is observed. However, this time-depen- 
dent flow does not include that which is rotating around the 
cylinder with a constant azimuthal velocity. This time-de- 
pendent convection can be buoyancy-driven, driven by the 
rotating magnetic field, or driven by some combination of the 
two forces. Indeed, as will be shown, there is a critical value of 
Ha for which Ra r = 0. 

Four distinct flow regimes have been identified. These four 
regimes are defined solely by the presence or absence of 
buoyancy-driven convection and time-dependent convection. 
No attempt was made to subdivide the time-dependent flow 
regime (both oscillatory and turbulent) and quantify the 
various flow modes that can exist there. 


4.1 

Ra <Ra , Ra'; Ha <3.8 

The flow in this regime can be described as being in the 
laminar, time-independent convective state. The rotating 
magnetic field drives an azimuthal flow component with the 
flow velocity depending on the strength of the magnetic field. 



Experimentally, all of the thermistors in the liquid exhibit 
a smooth, time-independent reading. Although the thermistors 
do not register any temperature changes, a secondary meridi- 
onal flow may exist, indeed, as has been discussed previously 
(Davidson and Hunt 1987), a variation in the azimuthal force 
as a function of z, which may occur in a cylinder of finite 
height, can cause such a secondary flow. Also, because Ra < 
Ra c , the vertical temperature gradient is not strong enough for 
any buoyancy-driven flow to exist. 

4.2 

Ra c < Ra < Ra'; Ha < 1 

In this flow regime, the rotating magnetic field is relatively 
weak but still drives a laminar, azimuthal flow component. Ra 
is larger than Ra c , and so a buoyancy-driven meridional flow 
also exists. Figure 5 shows the transition between the regime 
described in Sect. 4.1 and the present one. The experimental 
conditions used to generate the transition shown in Fig. 5 are 
essentially the same as those used to generate the transition 
shown in Fig. 3. The major difference is that the data of Fig. 5 
were obtained in the presence of a rotating magnetic field with 
B 0 = 0.72 mT. Just as in the case when Ha = 0, above Ra L a 
single, non-axisymmetric roll sets in. However, this single, 
non-axisymmetric meridional roll rotates around the cylinder, 
driven by the rotating magnetic field. The L3 — LI and L4 — L2 
thermistor signals have a near-sinusoidal response as a func- 
tion of time. The L3 — LI signal is shifted 90 with respect to the 
L4 — L2 signal, and the maximum temperature measured at L4 
occurs before the maximum temperature measured at L3. The 
angular velocity and rotational period can easily be calculated 
from data such as is shown in Fig. 5. The angular frequency of 
rotation of the single, non-axisymmetric roll increases with 
increasing magnetic field strength. It should also be noted that 
the critical Rayleigh number at which the transition observed 
in Fig. 5 occurs depends on the magnetic field strength, with 
Ra c increasing as Ha increases. 


Rayleigh number 

1.4 1.5 1.6 1.7 Xl0 4 



Time (s) 


Fig. 5. Temperature differences of thermistors situated directly across 
from each other in the test cell versus time. The corresponding 
Rayleigh number is shown at the top of the figure. The rotating 
magnetic field strength is 0.72 mT 



Time (s) 

Fig. 6. Temperature differences between the bottom and top of the 
ceil versus time. The decrease in T b — T a at approximately 2000 s 
indicates the occurrence of the flow transition at Ra 1 

The data in Fig. 6 are from the same experimental run as 
those shown in Fig. 5. The difference in temperature between 
the top and bottom of the cell is plotted as a function of time. It 
can be clearly seen when the transition at Ra c occurs. At ap- 
proximately 2000 seconds, T b —T a decreases slightly. It is at 
this point that buoyancy-driven meridional convection sets in 
and causes convective heat transport between the bottom and 
top of the test cell. After the transition, A T exhibits a sinusoidal 
characteristic. The frequency of this signal is the same as that 
observed for the thermistor temperature differences shown in 
Fig. 5. The thermocouples embedded in the copper plates at the 
top and bottom of the cell are off-center. Therefore, the ob- 
servance of the sinusoidal signal indicates that the thermal 
conductivity of the copper plates is finite and that there is 
a small, time-dependent, radial temperature gradient across 
them. 

The angular velocity of rotation was measured at several 
different magnetic fields, while keeping Ra constant. These 
data are plotted in Fig. 7. Even at the lowest magnetic fields 



Fig. 7. Angular velocity of the meridional roll versus rotating 
magnetic field strength for Ra= 1.7 x 10 4 . The solid circles are the 
experimental data. The dashed line is a fit to the data assuming that the 
angular velocity of the roll is proportional to B 2 0 






used, magnetically-driven azimuthal flow was observed. How- 
ever, at sufficiently high magnetic fields, Ra c becomes larger 
than Ra , and the buoyancy-driven meridional flow ceases. At 
this point, the flow transitions to the regime described in 
Sect. 4.1. Thus, for Ra= 1.7 x 10 4 , the experimental data shown 
in Fig. 7 essentially cover the entire magnetic field range of this 
flow regime. The dashed line in Fig. 7 is a fit to the data as- 
suming that the angular velocity is proportional to B As 
described by Davidson and Hunt (1987), such a dependence is 
consistent with laminar flow. It should be noted, however, that 
what is plotted in Fig. 7 is not the angular flow velocity, but 
rather the angular velocity of the azimuthal rotation of the 
meridional roll. Indeed, no r dependence was observed in the 
azimuthal rotation of the meridional roll, as would certainly be 
the case for the angular flow velocity. 

4.3 

Ra > Ra", Ra f ; all Ha 

In this regime, the flow is in the time-dependent state. This 
flow can consist of two or more oscillatory modes, or, with 
increasing Ra, turbulent flow. Because Ra>Ra‘\ Ha is by 
definition not large enough to suppress buoyancy-driven 
meridional convection. This flow regime can be reached from 
that described in Sect. 4.2 by slowly increasing Ra until Ra> 
Ra f . Above the transition, the near-sinusoidal thermistor 
signals, such as are shown in Fig. 3, gain an additional 
buoyancy-driven time-dependent component. The occurrence 
of the transition can be clearly observed by taking the deriv- 
ative of the signals with respect to time, which results in data 
plots similar to those shown in Fig. 4. 

4.4 

Ra ; < Ra < Ra"; Ha > 1 

In this flow regime, Ra c is always larger than Ra’. The 
combination of buoyancy and magnetically induced forces is 
strong enough to induce time-dependent convection, but Ha is 
large enough to suppress buoyancy-driven meridional convec- 
tion. The effect of increasing Ra " as the rotating magnetic field 
increases is similar to the situation descdribed by Chandrasek- 
har (1961), where the Ra c of fluid in a rotating cylinder in- 
creases according to the Taylor number. Here, the external 
rotating cylinder is replaced by a rotating magnetic field, but 
the effects are similar. As Ha increases, Ra c increases. In 
general, the time-dependence of the thermistor signals is driven 
by both electromagnetic and buoyancy forces. However, time- 
dependent thermistor signals were always observed for Ha> 3.8, 
even with Ra = 0. At Ra f , which depends on both Ha and Ra , 
the magnetic field drives the fluid to a high enough velocity 
that it transitions from the laminar to time-dependent state. 
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Summary 

A cylindrical test cell was constructed to study liquid metal 
flows in a rotating magnetic field. The fluid flow was inferred 
from temperature measurements of thermistors immersed in 
liquid gallium. The measured value of Ra 1 at zero magnetic 
field was in good agreement with recent theoretical and 
numerical results. Immediately above Ra c , the thermistor 
signals indicated a signal non-axisymmetric roll. Such a roll 
was consistent with previous theoretical predictions (Buell and 


Catton 1983) for a cylinder with an aspect ratio of 1, and was 
observed experimentally by Muller et al. (1984). With the 
magnetic field on, four distinct flow regimes were identified. In 
two of these regimes, laminar flow was observed. These are the 
flow regimes in which single crystal growth might be expected 
to be most useful. The results indicated that a rotating mag- 
netic field can increase Ra " and extend the range over which 
laminar flow is obtained. However, it should be noted that 
laminar flow could only be obtained for Hartmann numbers 
less than approximately 3.8. Above that value the flow was 
always observed to be time-dependent. A theoretical descrip- 
tion of the critical Rayleigh number's dependence on rotating 
magnetic field strength is the subject of future work. 
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